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Abstract 

We reexamine the production of J/d and other charmonium states for a variety of 
target-projectile choices at the SPS, in particular for the interesting comparison be¬ 
tween S-l-U at 200 GeV/c and Pb-|-Pb at 158 GeV/c as observed in the experiments 
NA38 and NA50 respectively. For this study we use a newly constructed cascade 
code LUCIFER II, which yields acceptable descriptions of both hard and soft pro¬ 
cesses, specifically Drell-Yan and meson production. This code divides the ion-ion 
collision into an initial phase involving hard interactions of the original nucleons and 
no soft energy loss, followed after the meson formation time by a ‘normal’ low energy 
cascade among the secondary particles. The modeling of the charmonium states dif¬ 
fers from that of earlier workers in its unified treatment of the hidden charm meson 
spectrum, which is introduced from the outset as a set of coupled states {ip, X*, V’^}- 
The result is a description of the NA38 and NA50 data in terms of a conventional, 
hadronic picture. The apparently anomalous suppression found in the most massive 
Pb-|-Pb system arises in the present simulation from three sources: destruction in 
the initial nucleon-nucleon cascade phase, use of coupled channels to exploit the 
larger breakup in the less bound x* ^|J' states, and comover interaction in the 
final low energy phase. 


1 Introduction 


The possible use of J/'0 suppression as a signal of unusual behaviour in rel¬ 
ativistic ion collisions, first suggested by Matsui and Satz [1], has attracted 
considerable experimental and theoretical study. Great interest has attached 
to the results obtained by the NASO collaboration for charmonium produc¬ 
tion in Pb-|-Pb collisions at 158 GeV/c: to the early findings presented at the 
Quark Matter 1996 meeting [2] as well as to the startling data later released at 
RHIG’97 [3]. The success of Glauber-like calculations of J/'0 production and 
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breakup in the p+A and S+U [4-6] systems, coupled with a failure of Glauber 
to provide an equally good description of the apparently accelerated absorp¬ 
tion in Pb-|-Pb has been widely interpreted [2,3,5] as a signal of QCD plasma 
creation in these collisions. The very sharp behaviour of the J/'^ yield as a 
function of transverse energy Et seen in the later experiment [3] has especially 
attracted attention. 

We attempt to retrace this ground theoretically, employing a new, two phase 
cascade approach, described in detail elsewhere [7,8], combined with a varia¬ 
tion of the Satz-Kharzeev model for production and annihilation of charmo- 
nium in the initial baryonic collisions. This modeling described below, allows 
the coupled-channel aspect of the hidden charm spectroscopy, {'?/’, X*, to 
play a more central role. The comparison with Glauber theory based models 
is done with two purposes in mind: first to understand the differences with the 
cascade if any, and second to help place the cascade on a firmer foundation, 
paradoxically, by indicating how similar the hrst high energy phase of the cas¬ 
cade is to the Glauber model. In this hrst application, we include partons in 
a minimal fashion, to describe for example Drell-Yan production. Hence we 
are testing a ‘purely’ hadronic description of the anomalous Pb-|-Pb measure¬ 
ments. From the evidence presented in hgures (10), (14), (15), it would appear 
such a test is justihed. 

It has been pointed out that a hadronic picture might succeed [6] without 
invoking quark-gluon plasma (QGP) creation, if at least part of the seemingly 
anomalous suppression in Pb-|-Pb could be produced by comover annihila¬ 
tion, i.e. by interactions of the J/'^ with secondary mesons generated in the 
ion-ion collision. The second phase in LUGIFER II, which is a low energy 
cascade, perforce includes the effect of J/'^ destruction through such comover 
rescattering. 

We begin with a description of the motivation behind LUGIFER II and a 
brief outline of the two step cascade. We attempt to separate hard and soft 
processes by time scale, so as to permit partonic and hadronic cascading to 
be joined naturally, in a modular fashion. The separation is effected through 
the use of a short time scale, automatically present at high energies: the time 
Tab taken for the two interacting nuclei A and B to traverse each other in 
the global collision frame. The uncertainty principle allows hard interactions 
involving sufficiently high energy-momentum transfer, i.e. for Q~^ < Tab, 
to take place in the hrst and very rapid cascade. Soft processes involving low 
tranverse momentum are not completed until later. Thus in the initial fast 
cascading the nucleons lose no energy but are still aware of the number and 
nature of the two-particle collisions they have undergone. 

Specihcally, the method [8] consists of running the cascade in two stages. The 
hrst is a high energy fast-time mode in which collision histories are recorded 
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and fast processes (Drell-Yan and charmonium production) are allowed to 
occur. Using the entire space-time and energy-momentum history of this stage, 
a reinitialisation of the cascade is performed using elementary hadron-hadron 
data as a strict guide. The hnal positions and momenta of baryons in the 
hrst phase, and the number of collisions they suffer are recorded and used to 
generate produced mesons together with their initial momentum and space- 
time coordinates. In the initial ion-ion collision the interacting nucleon paths 
are almost along light-cones. The second cascade begins at TaB) the time 
of the last nucleon-nucleon collision, with initial conditions specihed by the 
reinitialisation, but no secondary interactions are allowed until a formation 
time for produced mesons has passed. The participants in the second phase 
are generic mesons, thought of as of gg-like in character with masses centered 
near Mqq = 700 MeV and in the range Mqq ~ 0.3 — 1.0 GeV. Generic baryons 
consisting of qqq are also included and are excited to rather light masses, 
Mqqq ~ 0.94 — 2.0 GeV [8]. All the generic hadrons decay via sequential pion 
emission. Normal stable mesons and baryons are also present, and terminate 
the decay chains. 


Many cascades [9-15,7] have been constructed to consider relativistic heavy 
ion collisions. Since the eventual aim of experiments designed to study such 
collisions is the creation of a regime in which the quark-gluon structure of 
hadronic matter becomes evident, it will ultimately be necessary to include 
the partonic degrees of freedom in such cascades. However, since at SPS and 
even at RHIG energies it is by no means clear that all initial or subsequent 
hadron-hadron collisions occur with sufficient transverse momentum to free 
all partons [16], at least a part of the eventual simulation must deal with 
collisions both of the initially present baryons, nucleons in fact, and hnally of 
all the produced mesons as well. 


Kharzeev and Satz [5] employ a model based on hadronic Glauber theory 
describing production and breakup of the J/'^ in ion-ion collisions, to demon¬ 
strate that such a picture cannot account for the degree of suppression seen 
in Pb-|-Pb collisions at the SPS. Reasoning similarly, we can make a close 
comparison of our treatment with their work. The required initial production 
of a cc pair is handled within an effective hadronic formulation both in our 
work and in that of Kharzeev, et.al. . There are, naturally, specihc and im¬ 
portant differences between Glauber theory and a cascade model, and it is 
partly these differences which permit the so-called anomalous suppression in 
Pb-|-Pb to be explained within a purely hadronic framework. Additionally, the 
cascade provides a real calculation of Et as a function of centrality, i. e. impact 
parameter, with no adjustable parameters available for producing agreement 
with the experimental measurements. Equally important, the interaction of 
cc states with comoving secondary mesons is treated dynamically in the soft 
phase of our two step cascade. 
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The overall degree of suppression in Pb+Pb, insofar as it differs from earlier 
work [5,6], results from a combination of effects; these are baryonic, coupled 
channel and comover in kind, with substantial contributions arising from both 
phases of the cascade. There are potential unknown variables: the production 
and dynamic time evolution of each charmonium state, the breakup proba¬ 
bilities against both baryons and mesons, the density of secondary mesons. 
This last is to a large extent predicted by the cascade, which must agree with 
actual inclusive hnal state meson and baryon distributions. There also exist 
constraints on the basic charmonium variables. The production is in principle 
determined in elementary nucleon-nucleon collisions, the baryonic breakup in 
p-|-A collisions. The rjjTi breakup cross-sections are not known directly from 
any measurements. If the relative 'ip production in iViV and ttN systems may 
be taken as a guide here, the 'ipTi breakup cross-sections could be expected to 
be directly comparable. We have used a factor 2/3 to relate the charmonium- 
meson to charmonium-baryon breakup cross-sections, but also employed equal 
charmonium-meson as a test. 

The success obtained in describing the meson spectra has already been pre¬ 
sented in [8], and of course is relevant to the degree of charmonium destruction 
by comovers. In particular, the selection of a meson formation time is tied to 
this latter issue. The differences between S-l-U and Pb-|-Pb, which exhibit a 
considerable increase in the product Ax B, arise from both baryonic processes 
and meson production. 

The source code to LUCIFER II is made available on the world wide web at 
http://bnlnth.phy. bnl.go^ , and may be downloaded either directly from the 
web site or by anonymous ftp to bnlnth.phy.bnl.gov. It is in C and should be 
relatively easy to port to any UNIX system. Linux, AIX, SunOS, IRIX and 
HP/UX ports have at one time or another been made. 


2 The Two Phase Cascade 


We present in this work a mere outline of the cascade architecture, details hav¬ 
ing been provided in earlier work [8]. We already noted the global time scales 
which divide the cascade into two steps, the hrst loosely designated as ‘hard’, 
the second as ‘soft’. This separation has been discussed often in the literature 
[17-21]. Energy loss and meson production associated with low transverse mo¬ 
mentum pt are slow, ‘soft’ processes. In contrast stand fast or ‘hard’ processes, 
involving large pt, of which production of high mass Drell-Yan pairs [22-24] 
is a good example. The A-dependence of minimum bias Drell-Yan data [23] 
in p-|-A collisions, see Fig(l), suggests that high mass lepton pair production 
occurs only at the highest collision energy. For a theory of charmonium sup¬ 
pression to be taken seriously, it must simultaneously explain this striking 
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feature of Drell-Yan data, and the substantial soft energy loss experienced by 
the projectile nucleon. If not, an A-dependent suppression in p+A might be 
built into the model, attributable to production in successive projectile-target 
nucleon-nucleon collisions occurring at lower and lower energies, where the 
chance of J/'^ production is considerably less. The present model does well in 
this regard as seen in hgure 1. 

Nevertheless, calculations with a purely hadronic cascade [7] very well describe 
the energy loss and inclusive pion spectrum seen in massive Pb-|-Pb collisions 
at SPS energies, see Fig(2). These apparently contradictory features in fact 
can be united into a resonance based multi-scattering scheme. Following the 
high energy cascade stage in which collision histories are recorded and hard 
processes engage, the cascade is reinitialised and a second hadronic cascade 
is carried out at greatly reduced energy. Figure 3 shows the hnal positions 
of baryons in the first phase, and the almost light-like paths for the particles 
engaging in the initial cascade. 

In the reinitialisation, the nucleon-nucleon interaction history, together with 
the entire trajectory of each participant, is used to set up groups of nucle¬ 
ons which have mutually interacted in the hrst stage. The experimentally 
known averages, and the fluctuations inherent in NN scattering, in energy 
loss, multiplicity, and character (flavour, etc.) of produced mesons are all used 
to produce additional generic mesons associated with each nucleon group. 
Momentum, charge, baryon number and flavour conservation are imposed on 
each group of baryons and associated mesons. These generic mesons are the 
principal cascaders in the second step, along with the baryons. 


2.1 First Phase: High Energy Collisions 


As indicated, the procedure used here is relatively straightforward and in 
outline resembles the eikonal or Glauber calculations made by previous re¬ 
searchers [5,6,4], but retains the random, fluctuating, collisional nature of a 
cascade. This stage serves to establish the space-time geometry of the inter¬ 
actions between the target and projectile nucleons. Any actual hard processes 
which occur in this stage result in real energy loss, but soft processes do not 
occur yet, they are delayed until the reinitialisation. So, if a Drell-Yan pair is 
produced, its energy is immediately subtracted, similarly for a cc state or an 
open charm pair (we do not produce open charm at the moment, except by 
breakup of hidden charm mesons). Clearly jets or mini-jets (partons) could be 
produced in the same way and their hard evolution followed up to hadronisa- 
tion, when the products would simply be included during the reinitialisation 
before the second phase. 
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In practice, since Drell-Yan and J/'^ production are very rare processes, it is 
not necessary to actually globally conserve energy in making them, and this 
would in any case lead to considerable difficulty since the basic production 
must be considerably increased artihcially to make calculations possible in 
a hnite time. Therefore, energy is conserved in the hrst stage in the sense 
that realistic rapidity and pt (and mass for Drell-Yan pairs) distributions are 
employed for the appropriate elementary NN collision energy. 

Needless to say, the hrst stage interactions of charmonium preresonant states 
produced in the hrst stage cannot be neglected. This is because the cc pairs are 
produced relatively quickly, and may undergo high energy collisions with tar¬ 
get or projectile nucleons still in the way. Therefore these interactions are also 
counted for each charmonium, though actual breakup or feeding to another 
channel is left until the reinitialisation. 

The two phase approach has a highly benehcial ehect on frame dependence 
of the cascade, which has been touted by some authors as a possible limiting 
factor on the utility of cascade models for the description of ultra-relativistic 
ion collisions. Preliminary calculations at RHIC energy, ^/s = 200, for a central 
Au-|-Au collision, in our context a worst case scenario, yield at most a 10-15% 
effect on {dN/dy)^,- between the extreme cases of a cascade run in the global 
LAB frame and one in the global CM frame. Frame dependence at the SPS, to 
which energies the present work is limited, is negligible within statistical errors. 
The reason for the great reduction in frame dependence is simple and generic. 
The first, high energy cascade leads, see Fig(3), to particles travelling uniquely 
on light cone paths. Any residual frame dependence then enters only through 
the second phase cascading, which involves greatly reduced energies and is 
not very large. These results suggest that while frame dependence might be a 
limiting factor on pure hadronic cascades at LHC, it is certainly not a major 
difficulty at RHIC or SPS. But even partonic cascading is subject to similar 
problems, through the use of albeit reduced but still hnite cross-sections. 


2.2 Reinitialisation 


The fast cascade history is used to set up initial conditions for the second low 
energy cascade. To begin, one needs positions and momenta for the baryons 
and for the mesons expected to be produced from the initial baryon-baryon 
collisions. As we have noted already, the baryons are formed into groups which 
have collided with each other in the hrst stage. The structure of the groups 
is virtually dictated by consideration of the p-|-A system, where the projectile 
proton collides successively with those target nucleons which are directly in 
its path. 
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This simple grouping can be easily generalised to A+B collisions by using a 
procedure which kinematically symmetrises each group with respect to target 
and projectile. One first selects a nucleon which had the maximal number of 
collisions from all nucleons, both target and projectile, and adds to the group 
all those nucleons with which it collided. Since only high energy collisions are 
involved in the first stage, these colliders are all going (in an equal velocity 
frame) in a direction ‘opposite’ to the originally chosen particle. In p+A one 
is hnished at this point. For A+B one proceeds to choose from among the ‘op¬ 
posite’ colliders just added that nucleon which itself had the maximal number 
of collisions. All the nucleons with which it collided are hnally added to the 
group. Although geometry has carefully been ignored in this construction, it 
clearly plays its accustomed role. 

There is no doubt some limited freedom involved in constructing the groups. 
Quantitatively however, at least some of the possible alternative methods make 
little difference to the final results. This comes about not least of all because 
the hnal cascade phase really occurs at much lower energy, even in much higher 
energy RHIC collisions. The model resembles an extension of the wounded 
nucleon model [25]. We have essentially marked each incoming nucleon with its 
multifaceted history, and one might refer to this as a ‘painted’ or programmed 
nucleon model. 

The present work adds to the interaction choices the possibility of producing 
and destroying hidden charm cc states in each of the two steps of the overall 
cascade. The production is accomplished almost totally in the high energy 
phase, and uses elementary production cross-sections normalized to pp mea¬ 
surements. Breakup is done using the collision history for the charmonium 
collected in the hrst stage, through an interaction matrix described in the 
section on coupled channels. These processes are constrained by pp and p+A 
measurements [23,2,26]. 

The hnal step in the reinitialisation places mesons and baryons in position 
and time to restart the cascade. The four-momenta, in the global frame, of all 
particles are known, as are the space-time coordinates of the initial nucleons. 
It was thought best to distribute mesons produced within a group randomly 
along the paths of the baryons in that group. This choice would seem reason¬ 
ably consistent with locality principles. In any case, sins committed in this 
way are remediated by the formation time that must be attributed to each 
meson before it can begin to interact in the last phase. See Fig(3) for the hnal 
distribution of baryons following the hard cascade, and Fig(4) for a spacetime 
picture of the collisions in the entire cascade, both soft and hard phases. 
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2.3 Soft Phase 


The second phase then involves generic resonances, baryonic and mesonic, like 
the A, iV*, p and tt in quantum numbers. As stated, these have masses be¬ 
tween Mat = 0.939 GeV and about 2 GeV for baryons and from 0.3 GeV to 
1.1 GeV for non-strange mesons, and appropriately higher for strange mesons. 
These produced resonances are not allowed to interact until a sufficient forma¬ 
tion time has transpired. This time, rj, is a real adjustable parameter in the 
model, to be hxed perhaps from p-|-A, or as in earlier work from light nucleus 
collisions (S-l-S) [7]. We imagine that these are broad s-wave excitations of the 
underlying representatives, not series of special and very narrow states such 
as are tabulated in the particle data book. These narrow states cannot be ex¬ 
cited very much in actual ion-ion collisions, and do not carry much of the pp 
cross-section. Finally, these generic resonances decay by sequential pion emis¬ 
sion into lower mass excitations, losing mass with each decay until physical 
tt’s, p’s, and iF’s, or stable baryons are the only open channels. 

The decay time for the generic excitations is a second real parameter, per¬ 
haps to be taken inversely proportional to the excitation mass, but here for 
simplicity hxed at ~ 1/125 MeV“^. It is often assumed that one’s lack of 
knowledge of resonance-resonance scattering opens up to cascade models a 
deep well of adjustable parameters. This is not the case here. We employ, as 
did Gottfried [17], a ‘universality’ principle for soft interactions. Surely, for soft 
baryon-baryon interactions all hue details, excepting perhaps size and mass 
thresholds, are irrelevant in determining what must be, basically, interactions 
driven by many gluon exchange. We ignore size differences for the moment, 
with the exception later, of such very small objects as charmonium mesons. 
This limits the number of free parameters in the model to a minimum, in fact 
only a few, so far the two times Tf and r^. The important pp, pp and np data, 
determined from experimental measurements over a wide range energies at 
which cascading takes place, are the primary inputs, and serve to hx all of the 
baryon-meson and baryon-baryon interactions. Meson-meson interactions are 
hxed by appealing to constituent quark model counting. 


3 Model for Hadron-Hadron Interaction 


The objective of the cascade approach to ion-ion collision is to proceed from 
a knowledge of elementary hadron-hadron collision to a prediction of the far 
more complex many body event. Many approaches have been put forward 
[11,13,14] including strings [9,10], but we prefer to retain a particle nature for 
the cascade. 



The required input is a model for the elementary hadron-hadron system, be¬ 
ginning with nucleon-nucleon but easily extended to meson-nucleon and ulti¬ 
mately applied to any two body hadron-hadron collision. The basic processes 
are elastic scattering and inelastic production of mesons. The latter we divide 
into the well known categories [27]: diffractive scattering, referred to as single 
diffractive (SD), and non-single diffractive (NSD) [28]. A graphic description 
of these processes is given in Fig(5). The SD process, leading to a rapidity 
gap between one of the leading hadrons and the produced mesons cluster, is 
associated with the triple Pomeron coupling [29], while the NSD production is 
attributed to single (and double) Pomeron exchange and presumably results 
in the observed meson plateau. These diagrams are the basis for our hadron- 
hadron model but must be supplemented by an intermediate picture which 
allows us to apply them, not only to hadron-hadron interactions in free space 
but also inside a nuclear environment. 

The generic mesons depicted in Fig(5) and generic baryons, having rather light 
masses selected in the ranges suggested above, constitute the basic elements 
for rescattering in the second-phase cascade. Our principal phenomenological 
sources are tied to the rapidity [27] and multiplicity [28,30-32] information 
obtainable from the elementary collisions. Care must be taken to describe the 
observed multiplicity distributions. We have chosen to impose KNO scaling 
[7,33] in our parameterisation of these distributions, despite the relatively 
small deviations seen for example in the CERN UA5 experiments [28]. We use 
the distributions obtained in reference [7], covering a wide range of energies 
from a few GeV to lOO’s of GeV for pp and finally to 200 GeV/c and 1.8 
TeV/c for pp [28,32]. 


4 Coupled Channel Model for Charmonium 


The treatment of the hidden charm cc mesons within a ‘purely’ hadronic 
code presents some problems, perhaps not fully solvable within the effective 
hadronic treatment of such states. We do not deviate much in spirit from the 
work of previous researchers [4-6], but the devil lies sufficiently in the details 
to produce some quantitative effects. The production of charmonium mesons 
is almost completely limited to that coming from nucleon-nucleon collisions 
at the highest energies, i.e. in the initial high energy cascade, not by hat but 
by the greatly reduced collision energies in the second phase. Destruction of 
the charm meson precursors, in contrast, can occur in the hrst baryonic phase 
and also later in collisions with generic mesons and baryons in the second, 
low energy phase, i. e. on comovers. It is in the destruction of the charmonium 
states we differ most, ascribing a more direct role to the presence of the higher 
mass X and pj' mesons, for which in fact breakup is far easier. We include in 
Fig(6) a level diagram showing the relevant charmonium states to make the 
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picture as clear as possible. 


In the actual calculations of the above cited references both production and 
breakup are treated as instantaneous. There is no J/'^ formation time in 
the high energy phase or its near cousin, the Glauber or eikonal modeling. 
Kharzeev et al. [5] in fact justify such a choice by referring to microscopic 
production of charmed quark pairs, the subsequent formation of a preresonant- 
resonant state from which all charmonium mesons emanate, and the breakup 
in relatively hard scattering by gluons radiating from nearby nucleons. All of 
this can be incorporated in the initial cascade by ascribing to NN collisions 
a production cross-section normalized to the measured elementary pp and pn 
data and similarly a baryonic breakup deduced from say p-|-A measurements. 
Nevertheless, the eventual result is an effective hadronic modeling for the char¬ 
monium states, ascribing to a produced cc pair a hadronic state which can be 
destroyed. 

We imagine that the primordial cc pairs are originally produced essentially in 
plane wave states. Clearly, both singlet and octet color states are involved. 
This view would seem to be a reasonable one given the predominance of open 
charm production over hidden charm production in free space NN collisions. 
We further suppose that in elementary NN collisions the cc pair eventually co¬ 
alesces, with a state dependent probability, into a J/'^, or y. The time which 
elapses will be determined by the size of the bound state and the probability 
that a transition occurs. The probability of formation will depend critically on 
the relative momentum of the coalescing pair as well as on their spatial sep¬ 
aration. In any picture of charmonium generation there must needs be some 
formation time Th for the bound state, which may be, in general, longer than 
the total duration of the fast, baryonic cascade Tab- Then, in an ion-ion colli¬ 
sion, additional transitions may be induced: into the continuum, i.e. breakup, 
or perhaps also between bound states. 

Therefore, whether one sees the early evolution of the eventual charmonium 
as a preresonant state or as plane waves may be immaterial. Given the small 
size of the J/'^ as opposed to the much larger ip' and y, the separation of the 
c and c in the plane wave picture could equally well serve as a distinguishing 
feature. What can differentiate our calculation from the earlier models is the 
possibility of transitions between charmonium states dynamically in a nuclear 
environment. Gertainly, the y^, y^ and y^ states are produced considerably 
more copiously in basic pp collisions, with perhaps as high a ratio as y/'^ = 
4 — 5 [34], and they decay appreciably into J/t/’, with branching ratios in the 
range Tb/T ~ 12 — 25% [35]. The tjj' also feed some 57% into J/-^. It follows 
that one cannot ignore their presence. 

This point becomes even more signihcant when one considers what the breakup 
probabilities for the higher mass charmonium states are likely to be, either 
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in the fast or slow cascades. These heavier objects are considerably larger 
spatially and might well have total cross-sections on baryons or mesons pro¬ 
portional to the square of their colour dipole radius [36]. In any case, in the 
precursor states of the initial rapid cascade the spatially larger charmonia will 
have more time for collisions before forming, and should then possess larger 
effective breakup strengths. It is interesting that the J/'^ production cross- 
sections on baryons and on mesons differ little, perhaps also a reflection on 
annihilation. We interpret the observed suppression of J/'^ as a measurement 
of the tjj' and/or y strength as well. One of our conclusions will be that a 
considerable portion of the anomalous suppression seen in Pb-|-Pb, even for 
quite large impact parameters, is a result of breakup in the higher charmonium 
states and an extinction of their free space feeding down. 

What extra parameters has our model introduced relative to other treatments? 
We introduce a reaction matrix for charmonium states Rij which permits 
transitions between the states as well as diagonal, breakup, elements. Unitar- 
ity constrains this matrix. The diagonal elements should be present for other 
practitioners also but are not, in general, since only one preresonant charmo¬ 
nium state is usually considered. The only signihcant off-diagonal component 
for our purposes is that for the ip to tp' transition, shown below to be small, but 
potentially influential. The quantitative reaction matrix Rij is specihed in the 
ensuing sections treating the explicit calculations. Indeed, in our “standard” 
calculation, encompassing all major results presented below, this matrix is lim¬ 
ited to diagonal in form, thus reducing the number of adjustable parameters. 
The breakup matrix elements for the second phase cascade, i.e. for breakup 
on comoving mesons, are scaled from the first phase by counting constituent 
quarks. 

The formation time for secondary mesons, tj, controls the initiation of the 
second cascade and thus the onset of comover destruction of charmonium and 
the density of comovers. A reasonable choice for this parameter is r/ ~ 0.5 — 1 
fm/c, and this is in fact consistent with the production of tt mesons in ion-ion 
collisions at the SPS [8,37,38]. The high density of comovers which obtains at 
these times implies they play a considerable role. The effective formation time 
is actually somewhat longer, since it is increased by the duration of the fast 
cascade, i.e. Rff = tj + TabI‘^- 

The energy dependence of the elementary i\/ip production cross-sections is 
shown in Fig(7). The sharp dependence of aj/^ on energy near the SPS values 
y/s = 17 — 20 GeV implies that virtually all production occurs in the high 
energy phase. Drell-Yan exhibits a similar behaviour. 
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5 Drell-Yan 


The high energy phase, designed to record the initial interactions of the nucle¬ 
ons in the two colliding nuclei also provides the basis for our estimate of mas¬ 
sive dilepton production, i.e. Drell-Yan, an important side of the quandary 
we faced at the start. We limit ourselves to the canonical FNAL [23] p-|-A 
measurement at 800 GeV/c, but in fact the method of calculation guarantees 
agreement with the lower energy p-|-A and A-l-B collected by NASO [2]. The 
calculation is straightforward and was a first example, for us, of a hard process 
added to the initial high energy cascade. Drell-Yan is generally considered to be 
calculable perturbatively for dilepton pairs with masses in excess of = 4 
GeV. Production in the short time dehned by such masses proceeds without 
energy loss and leads to the A-dependence shown in Fig(l). The E772 results 
[23] imply very close to linear A dependence. To perform the Drell-Yan mi¬ 
croscopically we have introduced parton structure functions [7], and of course 
the requisite parton variables, into the code. But the curves in Fig(l) could 
have been obtained purely geometrically from the elementary production rates 
and the high energy phase only; very little production comes from the second, 
low energy phase. Any cascade which does not correctly describe this feature 
of Drell-Yan is in danger of producing spurious charmonium suppression by 
means of premature energy loss. Drell-Yan and J/'^ production both can occur 
in second, third and higher order collisions of initial nucleons as well as in the 
hrst collision, and if energy is lost by the nucleons immediately following each 
collision, the result will be A-dependent ‘suppression’, since the production of 
all charmonium states drops sharply with decreasing energy. 

Later we introduce for purposes of comparison a survival probability for J/'^, 
which has as its denominator, aside from a nucleon-nucleon normalisation, 
the Drell-Yan yield corresponding to the charmonium yield in the numerator. 
Thus, apart from overall normalisation the ratio a{ip)/a{DY) and the 3/if) 
survival probability Pg = N^{survived)/N^{produced) are one and the same 
thing. This is because elementary Drell-Yan production and J/'^ production 
are treated in exactly the same way in our cascade, and the number of high 
mass lepton pairs tracks the initial number of charmonia. 


6 Charmonium Suppression in Nuclear Collisions 

6.1 Minimum Bias: p+A and comparison to Glauber 

We begin with the suppression in p-|-A for which comovers play little role, 
the product pA of atomic numbers is small in comparison to say Pb-|-Pb. 
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Even here, however, the first stage high energy cascade does not suffice for an 
accurate description, some of the suppression on baryons occurs only in the 
second stage, as slow J/'^’s emerging from the interaction region are caught 
by nucleons, or interact at low energy in the target. The nucleon-nucleus data 
provides a necessary constraint on the basic parameters to be used in baryon- 
baryon production and breakup. This simplihed system also provides a fruitful 
ground for comparison between the eikonal approach and LUCIFER. 

To facilitate a comparison with the cascade we have made our own calculations 
with the Glauber formalism, relying on the formula: 


dS, 


Gl 


d% ABa{NN 


da{AB) 

d^b 


= j d^s dz dz' pa{s, z) psip — s, z) Ia{.s, z) /b( 6 — s, U), 


( 1 ) 

( 2 ) 


Ia{.s,z) 


OO 

-(R-1) j dZAPAis,ZA)cTabs) , 

2 


( 3 ) 


for the differential survival probability oi J/ip produced in p-|-A collisions, with 
all integrals, including the one over b, carried out numerically to obtain the 
total survival probability. Here (Tabs is the J/'^ breakup cross-section and is to 
be determined from p-|-A data. We follow Reference [5] in this development but 
employ a simpler, hard sphere, version of the nuclear density pA{b, z) for the 
purpose of comparing Glauber and LUCIFER 11. It was instructive to extend 
this comparison to A-|-A collisions to demonstrate that even Glauber does not 
reproduce the canonical power law, implied in the experimental descriptions 
[2,26,37] which always are compared with a straight line ht on a log-log plot, 
supposedly arising from purely baryonic breakup. These results are displayed 
in Fig(8) for the J/'^ without coupled channels to make an easier comparison 
between cascade and Glauber. The J/'^ absorption cross-section is taken so as 
to reproduce the p-|-A observations at 800 GeV/c [23], aabs ~ 7.0 mb with the 
hard sphere conhguration. But this value is equally successful for the lower 
SPS energies. 

A second comparison, again for minimum bias production of J/t/’, appears for 
p-|-A in Fig(9). In this we use the coupled channel modeling, whose details we 
elaborate further now. The relative production of the different charmonium 
states is taken so as to reproduce the pp data from the ISR [34] for the y to 
J/ip ratio, i.e. x/t/’ ~ 4.5, and for appropriate p' production [39]. A general 
ballpark for the measured y to J/'^ production ratio in pp is a factor of 4 — 5 
[34,39] and our final results are rather insensitive to a choice in this range, 
since a decrease could easily be compensated by a small transition matrix 
element between J/'^ and y. 
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The iIj' to J/'^ production in pp is taken near 0.33 so as to reproduce the 
hnal J/'^ contribution from the eventual decay of the ip' seen in a variety of 
experiments at a range of energies [2,39], i.e. implying an eventual ip' to 3 /ip 
ratio ~ 0.15. The free decay of the ip' into 3/ip is of course given the standard 
[35] value 0.57. 

The (and all other charmonium states) are assigned their correct 

masses, to properly include threshold effects in breakup. But the branching 
into J/ip is taken the same 0.18, as a weighted average over the measured 
electromagnetic values together with a very small hadronic component, < 0.5% 
[35]. Again, small variations ~ 1 — 2% in this branching ratio have little effect 
on suppression in any of the charmonium states and can be compensated for by 
commensurate changes in in, for example, overall cross-section normalisation. 

A hrst inference to be drawn from Fig(8) for A-l-A and Fig(9) for p-|-A is 
that the hrst high energy cascade produces pure 3/ip dynamics very much like 
that in Glauber. For the above choices a aabs ~ 7.0 mb leads to very nearly 
the same yield with the cascade as with Glauber. A second lesson, key to 
our development, is that the coupled channel model reproduces the Glauber 
result for 3/ip, using a smaller direct breakup cross-section, here taken as 
o'abs{J/'ip) = 5—6.0 mb, and including indirect destruction via the considerably 
larger Cabsix) = ^ <^abs{J/'P) for y and perhaps higher for ip'. The increased 
spatial sizes of the higher states strongly support the use of larger absorption 
cross-sections. A hnal observation for collisions extending io A x B values 
encompassing S-l-U and further to Pb-|-Pb, is that pure Glauber theory and the 
hrst stage nucleon cascade, both produce lines curving appreciably downwards 
on log-log plots, deviating from any power law. One gathers there is a little 
bit of ‘anomalous’ suppression even in a bare bones, no comover, theory. 

One should keep in mind that the trne charmoninm states in p-|-A collisions 
are prodnced mostly ontside the nncleus and we and other workers are, for this 
hrst cascade, considering charmoninm progenitors, perhaps preresonances or 
perhaps jnst comoving cc pairs at a certain separation. The ehect is however 
the same as using ehective charmoninm states, instantaneonsly prodnced, as 
we and previons researchers [4-6] are doing. This sitnation is altered when 
one turns to the second stage cascade which begins later, when all mesons, 
charmoninm and others, may have had time to precipitate. 


6.2 Suppression in A+B Collisions 


We have, in the last section, already presented results of simulations for 
the high energy cascade alone, and demonstrated the close similarity of this 
hrst stage to an eikonal approximation. The p-|-A results, compared to both 
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Glauber and experiment, were displayed only for minimum bias collisions. To 
complete the picture one must allow the soft cascade to go forward, especially 
for ion-ion collisions where the production of mesons becomes very signifi¬ 
cant. A key parameter in second cascade is the delay, tj or more properly 
Teff, afforded by the time scale involved in the ‘soft’ formation of mesons. We 
reemphasize that we use a more or less standard value r/ ~ 0.5 —1.0 fm/c, but 
are constrained by the production of mesons in the most massive system con¬ 
sidered here, Pb-|-Pb. This production was considered extensively in reference 
[8]. Perhaps a 10 — 20% uncertainty might attach to Tf. 

There are two sets of data to be considered: first, minimum bias J/'^ cross- 
sections as a function of the product A x B of nuclear atomic numbers, and 
second the ratio of J/ijj yield to Drell-Yan yield as a function of centrality, or 
more specihcally transverse energy Ef. 

Our results for minimum bias are displayed for the combined effect of both 
cascade phases in Fig(lO). The anomalous suppression in Pb-|-Pb is well re¬ 
produced by the totality of our two step, but otherwise conventional, hadronic 
dynamics. Part of the additional suppression in Pb-|-Pb relative to S-l-U al¬ 
ready arises from the high energy cascade, coming from the increased y and tf)' 
breakup in the more massive nuclear collision. But a considerable differential 
suppression arises from comovers, some 40% of the difference between S-l-U 
and Pb-|-Pb, the rest arising from the hard cascade. Thus our conventional, 
hadronic explanation of the interesting NA50 measurements is multifaceted, 
being rooted in both of the two cascade stages, in the coupling of charmonium 
channels, and in the presence of comovers. Part of the anomaly however, is 
perhaps illusory in view of the ‘curving down’ seen for large Ax B in Fig(8) 
and Fig(9). 

The calculated minimum bias ip' suppression is compared to data in Fig(ll). 
The strong drop occasioned by the large increase from p-|-W to S-l-U or Pb-|-Pb 
is clearly present in the theory. As is evident in this figure the ip' breakup 
strength inferred from p-|-A proves sufficient for both S-l-U and Pb-|-Pb. 

The breakup cross-sections in these simulations are 6.6, 20.0 and 25.2 mb 
for the tp, y* and ip' respectively. These represent absorption in charmonium- 
baryon collisions, and are reduced by the constituent quark factor 2/3 in p)- 
meson. Variation of these meson-meson cross-sections upwards to full equality 
with charmonium-baryon leads to ~ 0.5% change in the overall ^/p) suppres¬ 
sion for Pb-|-Pb (see Fig(12)). This small change is muted by two circum¬ 
stances, even in the soft cascade the suppression from charmonium-baryon 
collisions is still 1/3, and also the presence of a saturation effect in the degree 
of ff/p) suppression. Hence, one obtains a robustness in the predictions. 

This perhaps surprising non-linearity for the ff/p) interaction arises because in 
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our model this state is intrinsically tied up with the higher states. Introducing 
off diagonal elements Rij would produce a family of solutions. We have left well 
enough alone; the present few modeling parameters, now mostly determined 
independently of S+U or Pb+Pb data, surely having produced already an 
adequate description of observation. 

The survival probabilities for the J/'^ in S+U are 0.50 and 0.87 in the hard 
and soft cascades respectively. The same hgures for Pb+Pb are 0.42 and 0.775. 
These results are both obtained using the model parameters described above 
and an effective formation time Te// ~ 0.6 + 0.9/2.0 fm/c, that is close to 1 
fm/c. Comovers play a signihcant role, increasingly more so in Pb+Pb, but 
one must restate that baryons as well as mesons represent comovers in the 
soft cascade. The ratio of breakup on mesons vs that on baryons is ~ 2 : 1 for 
Pb+Pb and less for S+U. 


6.3 Centrality: Dependence on Transverse Energy 


Perhaps the most striking features of the NA50 [3] measurements are contained 
in their plot of J/'^ suppression vs Et. Unlike the existing Glauber calcula¬ 
tions of transverse energy the cascade provides a built in E^ scale, which does 
not necessarily agree exactly with the experimental determination. There is no 
theoretical parameter to adjust the simulated spectrum and the comparison of 
theory and experiment is a further test of the cascade. NA50 plots for J/'^ and 
Drell-Yan show Et from neutral mesons (x3 presumably in order to include 
all charge states) within the pseudo-rapidity range rj = 1.1-2.3. To establish a 
calibration from LUCIFER II we first refer to their earlier Pb+Pb results [37] 
using a more central rapidity range rj = 2.1-3.4, and including both electro¬ 
magnetic (neutral) and hadronic calorimeters to estimate Et. This comparison 
is shown in Fig(13) [37], and indicates that LUCIFER II, with standard pa¬ 
rameters [8], provides a reasonable representation of the measurements. The 
small discrepancy between cascade and experimental endpoints, some 10%, 
should be kept in mind when examining the NA50 charmonium data, but is 
unlikely to have a major effect on the results. 

Charmonium breakup in the second phase cascade occurs both on baryons 
and on mesons, generic and stable. Since meson numbers can be large, in par¬ 
ticular for Pb+Pb, their contribution to suppression may be appreciable. For 
collisions of J/'^ with generic mesons of mass greater than the charmonium 
binding, i.e. the energy interval 5E ~ 600 MeV from the charmonium mass to 
the DD continuum, we employ a properly exothermic cross-section. We use 
aabs oc kf/ki for the enhancement of the reaction cross-section. The enhance¬ 
ment due to the inverse R factor is not large numerically however, yielding a 
couple of percent at most in the excess suppression attributable to comovers. 
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Figures (14) and (15) display the results of simulations for the two massive 
ion-ion collisions. The magnitudes use the calculated survival probabilities, 
normalised by the pp or p + D experiments. The rather low Et value at which 
the measured 3/if) suppression becomes pronounced obtains equally well in the 
simulation, and the same low level of J/'^’s is reproduced. The results reinforce 
the perception already created by the comparison with the minimum bias data. 
The hadronic two-step cascade is capable of describing the charmonium yields, 
J/'^ and ip' as well. Again, the source of the suppression is multifold, perhaps 
half from comover interaction and half from the hard cascade, both strongly 
influenced by our treatment of the charmonium states. The beginning of strong 
suppression in J/close to peripheral collisions, Et ~ 50 GeV, is a reflection of 
the role the heavier charmonium states play, and occurs at or near the correct 
observational point. 


7 Conclusions 


It appears that a conventional hadronic explanation of the minimum bias and 
central J/'^ and ip' suppression in A-l-B collisions is possible. This is accom¬ 
plished here with a cascade, not specially tuned for the charmonium sector 
alone, but consistent with soft energy loss processes, both meson and proton 
spectra, and with a broad range of Drell-Yan data. There are parameters in the 
model, notably the meson formation time, the generic meson decay constant, 
and certainly the elements of the charmonium reaction matrix. However, we 
have not made use of all of this freedom in obtaining the main results. We 
have used a pure diagonal reaction matrix to describe the coupled charmo¬ 
nium channels, no mixing was needed further than that implied in NN. The 
ratio of charmonium states produced in the elementary nucleon interactions is 
naturally somewhat uncertain, but strong hints are available from experiment 
and we followed these indications. 

Many of the other parameters are also constrained by p-|-A data or by reason¬ 
able assumptions, for example the relative size of breakup probabilities must 
be tied to the relative sizes of the charmonium states. The formation time for 
produced mesons is also strongly tied to inclusive meson production. Some 
authors [5] suggest theoretically that the ^/ip total cross-section on mesons 
must be drastically smaller than the 4-6 mb we use for breakup, but this 
argument is disputed by others [36], and it seems unlikely. Our direct experi¬ 
mental knowledge of the total and partial cross-sections, including any energy 
dependence, for J/ip or other charmonium mesons on the lower mass mesons 
is of course very limited. 

Evidently, the initial phase of our model, or that of others, for charmonium 
creation or breakup is not strictly nor conventionally hadronic. The production 
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and destruction of states in this phase cannot be of standard fully formed 
charmonium states alone. But whether one imagines these initial forms to be 
preresonant or a somewhat excited cc may not change the effective dynamics. 
We have chosen to view the early life of the hidden charm mesons as akin to a c 
and c in relatively narrowly dehned continuum wave packets and characterised 
by their relative distance. It is crucial that the cc pair be in both colour 
singlet and octet states. Then, despite the lengthier formation time for the 
more massive states, in fact precisely because of this, these states are harder 
to form, suffering more collisions. The dependence of an effective cross-section 
on radius then follows. 

For an elementary pp collision the eventual coalescence of this pair into a 
proper charmonium state could be described in a fashion following the for¬ 
mation of deuterons or other clusters [40]. In the case of pp there would be 
independent formation probabilities for the various spectroscopic states, but 
the reaction matrix would presumably be diagonal, assuming no signihcant 
hnal state interactions on the charmonium. The coalescence within a nuclear 
environment would be very different than in free space, with many more in¬ 
teractions preventing bound state formation. The more general matrix we 
proposed, Rij, would have a simple meaning in an extension of this picture 
to ion-ion collisions, its elements describing what happens to these early cc 
states as they propagate through the nuclear environment, above and beyond 
their fate in pure pp. 

What then has been learned about excited, dense, nuclear matter from the 
reduction in J/'^’s? Our earlier calculations [8] for a broader range of processes, 
suggested that very high baryonic and mesonic energy densities were achieved 
in central Pb-|-Pb interactions, pb ~ 6-7/(fm)^ and pe ~ 4-5 GeV/(fm)^ 
respectively and that these densities persist for quite long times r ~ 3-5 
fm/c, in the CM frame. Thus appropriate conditions for possible ‘plasma’ 
creation exist in the most massive collision. This matter density has been 
sensed in the theoretical comover breakup of charmonium, both J/'^ and ip', 
more so for Pb-|-Pb than for S-|-U. But should our model stand the test of 
time, and it has explained a good portion of existing data at the SPS, then 
the case for a non-conventional explanation is hard to establish as yet. It is 
always of course possible that partons actually are playing a less passive role 
than portrayed by hadronic modeling, that especially high gluon densities are 
achieved in the initial phase. If so, the use of purely hadrons seems to mock 
up the partonic behaviour rather closely, and even the early charmonium and 
Drell-Yan dynamics (consider, e.g. p-fA) seem to, effectively at least, £t into 
this picture. In any case our simulations do not rule out the creation of some 
form of partonic matter in ion collisions at SPS energies. They only make the 
necessity thereof less compelling. A microscopic, non-hadronic treatment of 
the internal strusture of the charmonium states might alter the entire picture, 
and may be necessary at higher energies. 
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At RHIC energies and beyond, surely even the soft hadronic measurements 
will be contaminated by contributions originating in harder parton interac¬ 
tions, degraded say, by later hadronisation. Nevertheless, it may be hard to 
hnd partonic footprints in the soft spectra. The generic mesons and baryons 
may simply be stand ins for mini-jets [16], effectively and correctly incorpo¬ 
rating divided parton cross-sections and hadronisation within hadronic colli¬ 
sions, decays and meson production. If appreciable partonic plasma is quickly 
formed, in a state close to thermalisation, the loss of energy to the baryonic 
background might very well be detectable. We intend to pursue such a path 
to hnding ‘signals’ of collective parton dynamics. 

The dynamics of charmonium suppression would seem ab initio to be similar 
at RHIC, but the greatly reduced duration of the initial hard cascade at = 
200 GeV might play a spoiler role. The quasi-hadronic modeling employed in 
this phase by all practitioners may be inappropriate at such increased energy. 
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Drell-Yan A-Dependence (Lucifer II) 

p+A @ 800 GeV/c 
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Fig. 1. A-dependence of Drell-Yan at 800 GeV/c: E772 (FNAL) vs LUCIFER. Min¬ 
imum bias dimuon production as shown is calculated microscopically using NAS 
structure functions [22] but could as well have been obtained directly from consid¬ 
erations of the total collision number. The same argument applies to the high mass 
dimuon cross-section as a function of Ef. 
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Inclusive Distributions: Pb+Pb @ 158 GeV/c 


NA49 Data vs Lucifer 11 
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Fig. 2. The calculated Pb+Pb rapidity spectra [7] at 158 GeV/c for vr” and protons 
compared to measurements by NA49. The latter are for total h~ and h'^ — h~ 
respectively 
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Baryon Final Positions 
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Fig. 3. Positions of baryons in Pb+Pb after first phase: Fast cascade. The time and 
longitudinal envelope indicates where the purely baryonic system reaches after the 
two massive nuclei pass through each other at this SPS energy 200 GeV/c. Although 
this envelope represents the end of the fast cascade, some charmonium-baryon col¬ 
lisions may take place later, i.e. in the comover phase. 
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Soft Cascade 



Incoming Projectile Incoming Target 

Fig. 4. Time evolution of the ion-ion collision. The distribution in space-time of 
collisions and decays in hard and soft cascades is shown for the minimum bias 
Pb-|-Pb system. The initial paths of the incoming nuclei are close to light-like and 
result in the dense initial blob of hard collisions. The soft cascade occurs after the 
formation time of mesons has passed, as indicated roughly by the constant proper 
time surface. 
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String-Like Model for Hadron-Hadron Scattering 


Elastic 




Fig. 5. Shown are graphic representations of the elements of the model for the 
elementary hadron-hadron collision: elastic, single diffractive (SD) and non-single 
diffractive (NSD). The meson groups introduced in both SD (with a rapidity gap) 
and NSD have a string -like character but divided into our generic resonances. It is 
customary to associate SD with a three pomeron coupling and NSD with one, two 
or more pomeron exchange. 
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Production 
V’/V ~ 0.6 


XM ~ 4.5 


Wr=0-i8 


v~ 1 


DD continuum 3738 MeV 



Fig. 6. Charmonium spectroscopy including higher mass states which are signif¬ 
icantly produced in pp and which feed strongly to the itp. Electromagnetic and 
hadronic decays of y* (a weighted average) and tjj' are both included in the indi¬ 
cated branching ratios. The production ratios are suggested by direct measurement. 
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C(pp->J/\|/) 



Fig. 7. Production of J/V’ from pp as a function of energy. The vrp cross-section is 
also known, and in fact is very similar to that for pp, but rarely plays a role with 
production generally significant only at the highest energies. 






/A|/ Survival Probability in A+A 


Glauber vs Cascade: Hard Sphere Density 


,. 0 - J 



- Glauber 

• Cascade 
- Power Law 


^'^ 10 ° lo' id id id id id 

A*A 


Fig. 8. Comparison for A+A between Glauber and cascade, the latter in a purely 
i/tp mode and both calculations employ ahr = 7.3 mb. The deviation from a power 
law is apparent for large Ax A. A hard sphere form is used for the nuclear density. 
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/A|/ Survival Probability in p+A 



A 


Fig. 9. Again comparison of survival probabilities for Glauber and cascade, the 
latter appearing both in pure J/V’ and coupled channel mode. A variety of absorptive 
strengths are illustrated. The use of a hard sphere density distinguishes these results 
from later cascade simulations; in particular leading to small changes in the breakup 
cross-sections. 


30 







/A|/ Suppression in pA and AB 

LUCIFER CERN SPS DATA 



‘proj'^cible 


Fig. 10. The whole range of yields for frompp and p+D to Pb+Pb calculated in 
the cascade and compared to SPS measurements at various energies. The absolute 
theoretical values are obtained by normalisation to nucleon-nucleon. 
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Fig. 11. Comparison of experiment vs simulation for ijj'. The Pb+Pb data from 
NA50 was rescaled to 200 GeV/c by the collaboration. The S+U data is taken from 
NA38. The cascade calculations, again normalised to nucleon-nucleon, reproduce 
the observed behaviour for p+A and the sharp drop in the if:' to J/V' branching 
ratios for the massive nuclear collisions. 
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/A]/ Suppresion 


2/3 



Fig. 12. Variation oi 3/if) suppression with the charmonium-meson cross-sections. 
We use 2/3 as the ratio to charmonium-baryon for the calculations in the paper. 
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Fig. 13. Transverse energy distributions from LUCIFER compared to experiment. 
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Fig. 14. Comparison between the cascade and NA38 transverse energy dependence 
for J/'i/', for S+U. 
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Fig. 15. Comparison between the cascade and NA50 transverse energy dependence 
for J/'ip. There are no discontinuities, of course, in the LUCIFER yields, but the 
general shape is reproduced. In both this plot and that for central S+U, the Et 
scale is that obtained directly from the cascade, no scale change has been made. 
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